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Abstract 

Long-chain flavodoxins, ubiquitous electron shuttles containing flavin mononucleotide (FMN) as prosthetic group, play an 
important protective role against reactive oxygen species (ROS) in various microorganisms. Pseudomonas aeruginosa is an 
opportunistic pathogen which frequently has to face ROS toxicity in the environment as well as within the host. We 
identified a single ORF, hereafter referred to as fldP (for //avoafoxin from P. aeruginosa), displaying the highest similarity in 
length, sequence identity and predicted secondary structure with typical long-chain flavodoxins. The gene was cloned and 
expressed in Escherichia coli. The recombinant product (FldP) could bind FMN and exhibited flavodoxin activity in vitro. 
Expression of fldP in P. aeruginosa was induced by oxidative stress conditions through an OxyR-independent mechanism, 
and an fldP-null mutant accumulated higher intracellular ROS levels and exhibited decreased tolerance to H 2 0 2 toxicity 
compared to wild-type siblings. The mutant phenotype could be complemented by expression of a cyanobacterial 
flavodoxin. Overexpression of FldP in a muff-deficient P. aeruginosa strain decreased H 2 0 2 -induced cell death and the 
hypermutability caused by DNA oxidative damage. FldP contributed to the survival of P. aeruginosa within cultured 
mammalian macrophages and in infected Drosophila melanogaster, which led in turn to accelerated death of the flies. 
Interestingly, the fldP gene is present in some but not all P. aeruginosa strains, constituting a component of the P. 
aeruginosa accessory genome. It is located in a genomic island as part of a self-regulated polycistronic operon containing a 
suite of stress-associated genes. The collected results indicate that the fldP gene encodes a long-chain flavodoxin, which 
protects the cell from oxidative stress, thereby expanding the capabilities of P. aeruginosa to thrive in hostile environments. 
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Introduction 

Microorganisms living in aerobic environments are constantly 
exposed to the harmful effects of reactive oxygen species (ROS), 
including H 2 0 2 and the superoxide radical, which are generated as 
unavoidable by-products of oxygen utilization [1]. In addition, 
commensal and pathogenic bacteria have to face the host oxidative 
response, such as H 2 0 2 production from phagocytes [1,2]. Aerobic 
organisms have evolved multigenic responses to prevent and/or 
repair the cellular damage potentially inflicted by these toxic 
compounds. Whenever defenses are overcome by the amounts of 
ROS produced, cells are afflicted by a condition called oxidative 
stress [1,2]. Protective mechanisms deployed by stressed organisms 
include regulation of membrane permeability, antioxidant and 



repair systems, and replacement of ROS-sensitive targets by 
resistant isofunctional versions. In microorganisms as distantly 
related as enterobacteria and cyanobacteria, induction of the mobile 
electron shuttle flavodoxin (Fid) appears to be a common feature of 
the antioxidant response [3,4]. Fids contain flavin mononucleotide 
(FMN) as prosthetic group and are largely isofunctional with the 
ubiquitous electron carrier ferredoxin (Fd), exchanging reducing 
equivalents with a promiscuous lot of donors and acceptors. Fid 
induction is assumed to act as a backup for Fd, which harbors a 
ROS-sensitive iron-sulfur cluster as redox-active cofactor and whose 
levels are down-regulated under conditions of environmental stress 
or iron starvation [5,6]. Accordingly, Fid overexpression has been 
shown to confer augmented tolerance toward various sources of 
oxidative stress in organisms with very different lifestyles, such as 
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Author Summary 

Coping with toxic reactive oxygen species (ROS) generated 
as by-products of aerobic metabolism is a major challenge 
for 0 2 -thriving organisms, which deploy multilevel re- 
sponses to prevent ROS-triggered damage, including 
membrane modifications, induction of antioxidant and 
repair systems and/or replacement of ROS-sensitive targets 
by resistant isofunctional versions, among others. The 
opportunistic pathogen Pseudomonas aeruginosa is fre- 
quently exposed to ROS in the environment as well as 
within the host, and we describe herein a new response by 
which this microorganism can deal with oxidative stress. 
This pathway depends on a previously uncharacterized 
gene that we named fldP (for /favoofoxin from P. 
aeruginosa), which encodes a flavoprotein that belongs 
to the family of long-chain flavodoxins. FldP exhibited a 
protective role against ROS-dependent physiological and 
mutational damage, and contributed to the survival of P. 
aeruginosa during in vivo infection of flies as well as within 
mammalian macrophagic cells. Thus, fldP increases the 
adaptive repertoire of P. aeruginosa to face oxidative 
stress. 



Escherichia coli [7] , rhizobia [8] and plants [9] . Unlike Fds, which are 
present in all major kingdoms, Fids are restricted to various groups 
of prokaryotes and some oceanic algae [10]. From sequence 
alignments and structural considerations, they can be divided into 
two classes, short-chain and long-chain Fids, which differ by the 
presence of a 20-amino acid loop of a so far unknown function [1 1] . 
Phylogenetic analyses indicate that the two lineages have diverged 
only once [12]. 

Pseudomonas aeruginosa is a free-living bacterium commonly found 
in soil, water, moist locations, and most man-made environments 
throughout the Earth. P. aeruginosa has a wide metabolic versatility 
as the reflection of a large and flexible genome with a substantial 
number of genes, which facilitates its adaptability to thrive in 
different habitats, and allows a quick response to diverse environ- 
mental stimuli and challenges [13]. This remarkable versatility 
enables P. aeruginosa to infect damaged animal tissues or immuno- 
compromised individuals, where it constitutes an important 
opportunistic pathogen highly prevalent in nosocomial infections 
[14,15]. Indeed, this bacterium is of particular concern to patients 
with cystic fibrosis (CF) who are highly susceptible to P. aeruginosa 
and suffer severe and often fatal chronic airway infections [16]. 

The present research seeks to identify and characterize putative 
long-chain Fids in P. aeruginosa which could play a protective role 
under environmental stress conditions. Both P. aeruginosa and its 
relative Pseudomonas putida contain a single-copy gene encoding a 
short-chain Fid [17,18], annotated as mioC due to the homology of 
the product with the homonymous Fid from E. coli, but no long- 
chain Fid orthologs have been so far reported in pseudomonads. 
We identified a gene, PA14_22540 (hereafter referred to as JldP, 
for flavodoxin from P. aeruginosa), which displays low but 
significant sequence homology to the Anabaena and E. coli Fid 
genes (named isiB and JldA, respectively), and whose recombinant 
product was able to display Fid activity in vitro. Expression of the 
fldP gene in P. aeruginosa was induced by H 2 0 2 treatment via an 
OxyR-independent pathway, whereas its disruption increased 
H 2 0 2 -induced killing and accumulation of intracellular ROS. The 
mutant phenotype could be complemented by transformation with 
the isiB gene from Anabaena. Overexpression of JldP mitigated 
H 2 0 2 -induced cell death in a mjrfT-deficient P. aeruginosa strain as 
well as the hypermutability caused by DNA oxidative damage. 



The presence of a functional JldP gene contributed to P. aeruginosa 
survival in two model systems of infection: cultured mammalian 
macrophages and Drosophila melanogaster. Improved bacterial 
endurance in Drosophila resulted in higher death tolls of the 
infected flies. In line with its presumptive adaptive role, the JldP 
gene was found to be part of a self-regulated operon belonging to 
the P. aeruginosa accessory genome, a collection of strain-specific 
gene clusters which are acquired en bloc and expand the genomic 
repertoire to fit the needs for survival in adverse environments. 
Then, the collected results indicate that the JldP gene encodes a 
long-chain flavodoxin which is induced when P. aeruginosa is under 
oxidative stress to exert a protective role against the physiological 
and mutational damage caused by ROS. 

Results 

Search and structural analysis of a long-chain flavodoxin 
in the P. aeruginosa genome 

We performed an in silico survey of putative Fids in the genome of 
P. aeruginosa PA14 and compared the retrieved sequences with those 
reported for Anabaena (IsiB) and E. coli (FldA) long-chain Fids, whose 
protective roles against oxidative stress have been extensively 
documented [12]. To find orthologs, the IsiB and FldA sequences 
were compared to the P. aeruginosa PA 14 genome using the Domain 
Enhanced Lookup Time Accelerated Basic Local Alignment Search 
Tool (DELTA-BLAST), which is sensitive in detecting remote 
protein homologs [19]. Four unique open reading frames (ORFs) 
displaying both sequence homology and similar domain organiza- 
tion were retrieved. They were a putative oxidoreductase with a 
covalently-linked Fid-like domain (PA14_58560), the repressor 
binding protein WrbA (PA14_51990), a short-chain Fid similar to 
E. co/i MioC (PA14_19660), and an ORF (PA14_22540, tentatively 
referred to as FldP), which displayed the highest similarity in length 
and sequence with the long-chain Fids used as baits. Analysis of 
ORF PA1 4_22540 indicated that it would encode a 184-amino acid 
protein with a molecular mass of ~20 kDa, which is in the range of 
those typically observed for long-chain Fids (170-185 amino acids). 
Multiple-sequence alignment between FldA, IsiB and FldP showed 
that while FldA and IsiB display 47% identity and 67% similarity, 
FldP exhibits 23% identity with both flavodoxins, and 50% and 
4 1 % similarity with IsiB and FldA, respectively. 

By using the JPred3 software we next performed a prediction of 
the secondary structure of FldP and compared it with those of IsiB 
(Accession number P0A3E0) and FldA (Accession number 
P61949). According to this prediction, FldP displays a significant 
similarity, at the level of the secondary structures, with both FldA 
and IsiB (Figure 1). The three proteins fit the common scheme of 
five (3-sheets intercalated with five a-helices, as well as the loosely 
structured region dividing |35 into (35a and (35b (between residues 
130 and 150 of FldP), which is typical of the long-chain class of 
Fids [11,12]. 

The fldP gene encodes a functional flavodoxin 

To determine if the product of the JldP gene is a functional Fid, 
the coding sequence was expressed in E. coli under the control of 
an inducible promoter (see Materials and Methods). The 
recombinant protein accumulated to high levels in the bacterial 
host, but unlike IsiB, which was readily soluble and assembled its 
prosthetic group in the E. coli cytosol [20], the P. aeruginosa protein 
was recovered largely as insoluble inclusion bodies (Figure SI A). 
Only after the simultaneous expression of a suite of molecular 
chaperones, a minor but significant amount of the protein could be 
solubilized and purified by affinity chromatography on Ni-NTA 
columns (Figure SIB). 
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Figure 1. Primary and secondary structures of FldP from P. aeruginosa. FldP sequence was compared with those of well-characterized 
flavodoxins from Anabaena (IsiB) and E. coli (FldA). Secondary structure predictions are represented above the sequence alignments. "E" identifies 
amino acids involved in p-sheets (yellow arrows), and "H", those which are part of a-helices (pink helices). The green arrow corresponds to the extra 
loop that characterizes long-chain Fids. Amino acids in bold (corresponding to W58 and Y95 in Anabaena IsiB) indicate the aromatic residues stacking 
onto the isoalloxazine ring system in IsiB and FldA, and their equivalents in FldP. 
doi:1 0.1 371 /journal.pgen.1 0041 63.g001 
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The purified protein showed a typical flavoprotein spectrum 
with absorption maxima at 374 and 446 nm, close to those of free 
FMN in aqueous solution (Figure 2A). The 446-nm peak, which 
corresponds to transition I of the flavin, is strongly red-shifted in 
the flavodoxins from Anabaena and E. coli (Figure 2A), indicating 
that the environment of the prosthetic group is more hydrophobic 
and/or less solvent-exposed in these flavoproteins. The flavin 
moiety of IsiB is sandwiched between the aromatic side-chains of 
Trp58 and Tyr95, in a coplanar conformation, and the Tt-7i 
interaction thus established has been regarded as the major cause 
for the spectral red shift [12,21]. The two amino acids are 
conserved in FldA, but not in FldP, where the tryptophan is 
replaced by a tyrosine and the tyrosine by a leucine (Figure 1). The 
substitutions will prevent aromatic stacking of the isoalloxazine 
ring system, and introduce conformational changes in its 
immediate environment, since both residues are located in flexible 
loops between 3(3-3ot and 4|3-4a (Figure 1). Then, the absence of 
red-shifted peaks in the visible absorption spectrum of FldP likely 
results from a combination of higher solvent exposure and 
decreased aromatic stacking on the flavin. 

FldP was assayed in vitro as a substrate of cyanobacterial 
ferredoxin-NADP(Fl) reductase (FNR). Figure 2B shows that 
purified FldP was able to mediate FNR-driven cytochrome c 
reduction in a concentration-dependent manner, with an apparent 
K M of 1.3±0.2 uM and a k, M of 22.0± 1.8 min" 1 . Under similar 
conditions, IsiB displayed a k vilt of 48.5 ±3.8 min 1 (data not 
shown). The collected results indicate that the product of the fldP 



gene displayed the structural and functional properties of a bona 
fide flavodoxin. 

Inactivation of fldP decreases survival of P. aeruginosa 
and increases intracellular ROS levels upon exposure to 
H 2 0 2 

In order to determine the functional role of FldP in P. aeruginosa, 
we tested the tolerance exhibited by a ^ZaP-deficient mutant strain 
to H 2 0 2 toxicity. In the absence of stress, the viabilities of the wt 
and fldP mutant strains were similar (Figure S2), but fldP cells were 
~5-fold more sensitive than their wt siblings to the H 2 0 2 
treatment [P— 0.024, Figure 3A). Complementation of this mutant 
with the fldP gene cloned in plasmid p2 (p2-fldP) led to an 
approximately 1.4-fold increase in the percentage of surviving cells 
respect to that observed with the wt strain, although the difference 
was not statistically significant (P— 0.3929, Figure 3A). Notewor- 
thy, expression of IsiB from the p2-uiB plasmid provided even 
higher levels of protection against the deleterious effects of H 2 0 2 , 
with a rise of 2.2-fold relative to the wt strain (P= 0.125, 
Figure 3A). 

Taking into account these observations, we further tested 
whether FldP may be involved in the control of intracellular ROS 
build-up, presumably the cause of H 2 0 2 -induced killing. ROS 
accumulation was quantified in bacterial extracts from the wt, the 
fldP and the complemented mutant strains after exposure to H 2 0 2 . 
Cells from the wt strain displayed a small increase of their ROS 
(-OOH) levels as the H 2 0 2 concentration was raised (Figure 3B). 
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Figure 2. FldP from P. aeruginosa is a functional flavodoxin. The 

FldP holoprotein was expressed in E. coli and purified by affinity 
chromatography as described in Materials and Methods. (A) Visible 
absorption spectra of purified FldP (black), FldA (blue), IsiB (green) and 
FMN (red) in 50 mM Tris-HCI pH 8.5. (B) Kinetics of cytochrome c 
reductase activity as determined by the increase in absorbance at 
550 nm. Numerals on the side of each curve indicate FldP concentration 
(in uM) in the reaction medium. The Fid-driven reaction sequence is 
indicated at the bottom. Experimental conditions are given in the text. 
doi:1 0.1 371 /journal.pgen.1 0041 63.g002 



Although this increase was not statistically significant, it was 
consistently observed in a number of experiments. On the other 
hand, lack of a functional FldP led to ROS build-up in the JldP 
mutant, whereas expression of FldP from a plasmid in the 
complemented cells decreased the total -OOH levels to those 
observed in untreated wt bacteria (Figure 3B). 

ROS accumulation was also detected in whole P. aeruginosa cells 
by using the fluorogenic dye 2',7'-dichlorofmorescein diacetate 



(DCFDA), and visualized by confocal microscopy. Figure 3C 
shows that the fraction of labelled cells above the detection 
threshold was significantly higher in the JldP mutant compared to 
their wt siblings (33% vs. 14%), but decreased to 2% after 
complementation with the rfi-jldP plasmid, presumably due to the 
effect of increased genie doses provided by the plasmid. 

Overexpression of fldP and isiB increases resistance to 
H 2 0 2 in a muf7"-deficient strain of P. aeruginosa 

It has been recently reported that P. aeruginosa strains deficient in 
the 8-oxodeoxiguanine system (GO) are particularly vulnerable to 
oxidative stress [22,23]. Specifically, mwiT-deficient cells showed to 
be the most susceptible to oxidants such as H 2 0 2 or methyl 
viologen. We therefore used this strain to further characterize the 
protective activity of FldP against ROS. A mutT P. aeruginosa strain 
was transformed with either nl-jldP or p2-isiB, and the resulting 
transformants were tested for their susceptibility to H 2 0 2 . Parallel 
controls were carried out by using the wt and mutT strains 
harboring an empty p2 plasmid. 

As shown in Figure 4, the mutT strain showed a 20-fold decrease 
in survival after exposure to H 2 0 2 , being significantly more 
susceptible than the wt strain (P= 0.01 19). However, when JldP or 
isiB were overexpressed in the mutT mutant, cell survival increased 
dramatically (8- and 16-fold, respectively, P— 0.01 19), relative to 
mutT transformed with the empty p2 vector (Figure 4). The results 
suggest that the antioxidant role of FldP and IsiB can partially 
compensate the increased susceptibility of ma/T-deficient cells to 
oxidative stress. 

FldP and IsiB decrease the mutation frequency induced 
by H 2 0 2 in a mut7"-deficient strain of P. aeruginosa 

DNA damage produced by ROS can lead to increased mutation 
frequencies [24] and emergence of adaptive phenotypes [25,26]. 
Particularly in mzrfT-deficient mutants, in which the damage 
produced by oxidative stress cannot be avoided, mutation 
frequencies can increase 100- to 1000-fold, leading to hypermu- 
tator phenotypes [22]. We therefore investigated whether JldP and 
isiB could display an antimutator effect and alleviate the 
hypermutability that is typically observed in a mjrfT-deficient 
background. Thus, we exposed a mu/Tstrain, which overexpressed 
JldP or isiB to H 2 0 2 and measured the mutation frequency by 
determining the emergence of mutants resistant to streptomycin. 
The mutT strain showed a ~ 1200-fold increase in the H 2 0 2 - 
induced mutation frequency (3.44x10 6 ), relative to the wt 
(2.83x10 9 ). Expression of either JldP or isiB in the mutT strain 
diminished the mutation frequencies to 4.49x10 7 and 
2.63 xlO" 7 , which represent a 13% (P= 0.05) and 8% 
(P= 0.0286), respectively, of the mutation frequency showed by 
mutT cells transformed with p2 (Figure 5). This last result indicates 
that overexpression of a functional Fid could avoid the majority 
(~90%), but not all H 2 0 2 -induced lesions produced as a 
consequence of mutT deficiency. The protective effect presumably 
results from the antioxidant properties of these flavoproteins. 
Accordingly, the antimutator effect conferred by both JldP and isiB 
was not observed when the spontaneous mutation frequency was 
tested (data not shown). 

Expression of the fldP gene is induced by H 2 0 2 treatment 
via an OxyR-independent pathway 

To gain further insight into the role played by FldP in oxidative 
stress tolerance, we studied the transcriptional induction of the JldP 
gene by H 2 0 2 treatment. The expression of JldP was monitored by 
semi-quantitative, two-step, reverse transcription-PC R (RT-PCR), 
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Figure 3. Role of FldP in cell survival and ROS accumulation upon exposure to H 2 0 2 . (A) Flavodoxin increases oxidative stress tolerance of 
P. aeruginosa exposed to H 2 0 2 . Cells from wt P. aeruginosa PA14 transformed with the empty plasmid p2, and from the isogenic fldP mutant strain 
harboring vectors p2, p2-fldP or p2-isiB were normalized to ~10 8 cells and then incubated with 50 mM H 2 0 2 for 30 min at 37'C. Viability was 
estimated as described in Materials and Methods. (B) Expression of FldP prevents ROS build-up in stressed P. aeruginosa cells. Bacteria were exposed 
to the indicated concentrations of H 2 0 2 for 30 min, and total peroxides (-OOH) were determined in cleared extracts using the FOX II assay as 
described in the experimental section. Measurements in (A) and (B) were carried out in triplicate for two independent experiments, and the results 
are expressed as means with their SEM. Statistically significant differences P<0.005 and P<0.05 are identified by ** and *, respectively (one-tailed 
Mann-Whitney test). (C) Fluorescent detection of ROS formation in cells from wt, mutant and fWP-complemented strains challenged with 25 mM 
H 2 0 2 for 15 min, using a ROS-dependent fluorescent probe (DCFDA). Panels on the left show membrane fluorescence after staining with FM4-64; 
those on the middle depict ROS-dependent fluorescence, with the merges displayed in the right panels. Bars show percentages of ROS-stained cells 
relative to membrane-stained cells. Values presented are the means ± SEM of 8-10 independent fields. Statistically significant differences respect to 
wt (P<0.0005) are identified by *** (one-tailed Mann-Whitney test). 
doi:1 0.1 371 /journal.pgen.1 0041 63.g003 
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Figure 4. FldP mitigates H 2 0 2 -induced cell death in a mt/f7"- 
deficient P. aeruginosa. The mufT mutant strain (transformed with the 
empty plasmid p2) was tested for resistance to H 2 0 2 , and compared 
with its wt parental strain (harboring p2), and with mutT mutants 
transformed with p2-fldP or p2-;s/B. Bacterial cultures were normalized 
to ~ 10 s cells and then incubated with 50 mM H 2 0 2 for 30 min at 37°C. 
Viability was estimated as described in Materials and Methods. 
Measurements were carried out in triplicate for two independent 
experiments, and the results were expressed as means with their ± 
SEM. Statistically significant differences (P<0.05) are identified by * 
(one-tailed Mann-Whitney test). 
doi:1 0.1 371 /journal.pgen.1 0041 63.g004 



Figure 5. FldP decreases the H 2 0 2 -induced mutation frequency 
in /77t/f7"-deficient P. aeruginosa. The muff mutant strain transformed 
with the empty vector p2 was tested for its frequency of mutations 
induced by H 2 0 2 (50 mM, 30 min), and compared with mutT mutants 
expressing fldP or isiB. The H 2 0 2 -induced mutation frequency was 
calculated as the ratio between the number of streptomycin-resistant 
cells and the number of viable cells and normalized to 1 in the mutT 
strain. The experiments were carried out in quintuplicate for three 
independent replicas, and the results were expressed as means ± SEM. 
Statistically significant differences (P<0.05) are identified by * (one- 
tailed Mann-Whitney test). 
doi:1 0.1 371/journal.pgen.1 0041 63.g005 



using expression of the constitutive housekeeping rpoD gene as 
control for equal amounts of cDNA in each reaction. Figure 6 
shows that expression of JldP in the wt strain was clearly induced 
(3.9 ± 0.7-fold relative to untreated cultures, P= 0.004), after H 2 O z 
treatment. 

We further investigated whether this increased expression of JldP 
was dependent on OxyR, the main oxidative stress response 
regulator of P. aeruginosa [2 7] . An oxyR deletion mutant strain of P. 
aeruginosa was constructed and tested for jldP expression in response 
to H 2 O z . Interestingly, no differences were observed in the oxyR 
strain compared to its isogenic wt, showing a 3. 0±0. 7-fold 
induction in the expression of JldP relative to untreated cultures 
(P— 0.0045, Figure 6). Importantly, we also evaluated expression 
of JldP in a second oxyR strain (ID54029) from the PA14 insertion 
mutant library [28], which yielded equivalent results (data not 
shown). These observations indicate that although JldP is a stress- 
responsive gene in P. aeruginosa, as in other bacterial species [3,4], 
this response is triggered by an OxyR-independent mechanism. 

FldP enhances P. aeruginosa survival within mammalian 
macrophages and during in vivo infection of Drosophila 
melanogaster 

The use of ROS to kill bacterial pathogens, such as H 2 0 2 
production from phagocytes, is a common feature of the innate 
immune response of eukaryotic organisms. Considering that FldP 
sheltered P. aeruginosa from oxidative stress under in vitro conditions, 
we tried to determine if this protective role of FldP could also 
provide an advantage to cope with the host immune defenses. 

As a first approach, we evaluated the capacity of the different P. 
aeruginosa strains to survive in the intracellular milieu of phagocytes 



by using monolayers of the macrophagic cell line RAW 264.7, 
which were inoculated with wt P. aeruginosa or its isogenic JldP- 
deficient strain, complemented or not with \>2-jldP. Then, with the 
addition of antibiotics to kill extracellular bacteria we were able to 
compare the proportion of bacterial cells of each strain that 
survived during a 3-h period inside the phagocytes (see Materials 
and Methods) by lysing the cell monolayer and plating the lysates 
on LB agar. Figure 7A shows that inactivation oiJldP produced a 
moderate but significant decrease of ~24% in the intracellular 
survival of P. aeruginosa in phagocytic cells (P= 0.0103). Impor- 
tandy, this decrease could be reverted by complementation with 
p2-JIdP, even surpassing the wt values. This result indicates that 
FldP is contributing to the intracellular survival of P. aeruginosa in 
macrophagic cells, probably by enhancing bacterial resistance to 
the ROS produced by phagocytes. 

To evaluate if this protective role of FldP could also be observed 
in the context of a whole organism infection, which is a more 
complex system than cultured cells, we used the insect-infection 
model of D. melanogaster, as previously validated to evaluate 
virulence traits of P. aeruginosa [29,30]. Thus, flies were fed with the 
same strains of P. aeruginosa mentioned above and left for periods of 
1 , 3 or 5 days, after which the ability of each P. aeruginosa strain to 
survive within the host was scored by counting colony forming 
units (GFU). Interestingly, the capacity of the JldP mutant to 
survive within the host dropped to 28 to 36% (P<0.05) at the three 
time-periods assayed, which could recover to wt values after 
complementation with p2-JldP (Figure 7B). 

These results suggest that increased bacterial loads could have 
an impact on host survival. To assess this, D. melanogaster flies were 
starved for 3 h and then continuously fed with P. aeruginosa. The 
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Figure 6. Induction and regulation of fldP expression. Total RNA 
was extracted from cultures (OD 600 = 0.8) of wt, oxyR and ID53714 
(PA14_22550 mutant) strains of P. aeruginosa which had been treated 
with 50 mM H 2 0 2 for 15 min or left untreated (NT). Semi-quantitative 
RT-PCR was performed with specific primers for fldP and the 
amplification products resolved by agarose gel electrophoresis. The 
rpoD transcripts were used as housekeeping controls. The upper panel 
shows relative expression levels averaged over at least 4 experiments 
(means ± SEM), with NT expression of the wt strain normalized to 1. 
Statistically significant differences at P<0.005 and P<0.05 are identified 
by ** and * respectively (two-tailed Student's f test). The lower panels 
illustrate typical results obtained with fldP and rpoD. 
doi:1 0.1 371 /journal.pgen.1 0041 63.g006 

number of surviving flies was monitored every day until all flies 
died. The results showed that the flies that had been fed with the 
wt strain of P. aeruginosa died faster than those fed with the JldP 
mutant (Figure 7C), which is in agreement with the better capacity 
of the wt P. aeruginosa strain to survive within the host. 
Complementation of the mutant bacteria with ^2-JldP increased 
the mortality rate to wt values (Figure 7C). 

Taken together, these findings suggest that FldP could be 
playing a role in P. aeruginosa pathogenesis by increasing its 
resistance to the ROS-dependent clearance carried out by the host 
immune system, which in turn would result in a more lethal 
infection due to a higher load of viable bacterial cells in the host. 

The fldP gene is a component of the P. aeruginosa 
accessory genome 

The P. aeruginosa genome is made up by a mosaic of "core" and 
"accessory and variable" gene clusters [31,32]. While the gene 
composition of the core genome is conserved in almost every strain 
of P. aeruginosa, genes belonging to the accessory genome are found 
in discrete patches, referred to as regions of genome plasticity 
(RGP), which can vary in occurrence and location among strains 
[33]. By using genome sequence information available online, we 
evaluated the occurrence of the JldP gene in thirteen P. aeruginosa 
strains (namely PAOl, PAH, 2192, 39016, C3719, LESB58, 
PACS2, M18, NCGM2.S1, B136-33, RP73, DK2 and PA7), whose 
genomes have been completely sequenced and made available at 
the Pseudomonas Genome Database (www.pseudomonas.com) [34]. 



The genome of the environmental strain Hex IT [35], which has 
been recently sequenced (Feliziani et al, unpublished), was also 
included in the analysis. The survey showed that only six of the 
fourteen strains (PAH, 39016, NCGM2.S1, B136-33, PA7 and 
HexlT) contained the JldP gene, suggesting that it is a component of 
the P. aeruginosa accessory genome. Indeed, JldP is located in the 
region of genome plasticity RGP32 [33]. With the exception of the 
taxonomic outlier PA7, all other strains showed a conserved synteny 
of RGP32, with/HP (PA14_22540) being the fifth gene in a six-gene 
cluster (PA14_22500 to PA14_22550 in the PA14 genome) 
(Figure 8). The DNA sequence of RGP32 is highly conserved 
among strains PA14, 39016, NCGM2.S1, B136-33 and HexlT, 
displaying up to 97% identity. The five ORFs accompanying JldP in 
RGP32 have been assigned putative functions, based on similarity 
with known genes (Table SI). Interestingly, RGP32 is flanked by 
two palindromic sequences (inverted sequence repeats) which could 
have played a role in the acquisition of this gene cluster. On the 
other hand, the genome of PA7 only showed orthologs for JldP 
(PSPA7_2449) with 62% identity, and for the two flanking genes of 
RGP32 (PSPA7_2618 and PSPA7_2620), both with 61% identity 
respect to their PA14 orthologs. Importandy, none of these genes 
share the genomic location observed for RGP32 in the other five 
strains. 

To further investigate the prevalence of JldP, we carried out 
PCR analysis using conserved primers (Table S2) on a collection of 
clinical and environmental isolates of P. aeruginosa, which were 
selected on the basis of being clonally different [36] . While highly 
divergent JldP versions could have been overlooked by this 
experimental approach, the high degree of sequence conservation 
of JldP observed among the characterized P. aeruginosa genomes 
makes this possibility rather improbable. The analysis showed that 
~23% of clones amplified a fragment corresponding to the JldP 
gene (Figure S3). Considering this low prevalence, the origin of 
RGP32 in P. aeruginosa is better explained as a consequence of 
DNA acquisition through one or more horizontal gene transfer 
processes rather than a genetic loss in those isolates which do not 
harbor RGP32. 

fldP is the first gene of an operon regulated by a LysR- 
type transcriptional regulator 

We subsequently analyzed the transcriptional organization of 
RGP32 in order to elucidate whether the JldP gene behaved as a 
single transcriptional unit or showed co-expression with other 
genes of the same variable region, thereby constituting an operon. 
We carried out PCRs using cDNA as template and primers 
designed to amplify fragments containing regions of two neigh- 
boring genes. We measured co-expression of those RGP32 genes 
which shared the same orientation and were located downstream 
from JldP JldP to PA14_22500). As shown in Figure 9A, we were 
able to amplify fragments between JldP and PA14_22530, 
PA14_22530 and PA14_22520, and PA14_22520 and 
PA14_22510. In contrast, no PCR fragments could be amplified 
between PA14_22510 and PA14_22500. Importantly, amplicons 
were obtained among all adjacent genes, even between 
PA14_22510 and PA14_22500, when genomic DNA templates 
were used as positive controls (Figure 9A). These results depict a 
transcriptional organization of RGP32 structured in one polycis- 
tronic operon containing JldP, PA14_22530, PA14_22520 and 
PA14_22510, and two monocistronic transcriptional units for 
genes PA14_22500 and PA14_22550. To strengthen further these 
observations, we performed RT-PCRs to measure transcripts of 
PA14_22530 and PA14_22500 in cultures which were treated or 
not with H 2 0 2 . Interestingly, expression of PA14_22530 exhibited 
a 3-fold H 2 0 2 -dependent induction, similar to that previously 
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Figure 7. FldP enhances P. aeruginosa survival within mamma- 
lian macrophages and during in vivo infection of Drosophila 
melanogaster. Comparison between the P. aeruginosa wt strain and its 
isogenic fldP mutant complemented or not with a functional copy of 
the fldP gene (p2-fldP), with respect to: (A) the fraction of bacteria that 
survived after 3 h in the intracellular milieu of RAW 264.7 macrophagic 
cells. Three wells of cells were used for each strain in each experiment, 
and all experiments were repeated three times. Statistically significant 
differences relative to the wt strain (P<0.05) are identified by * (one- 
tailed Student's f test); (B) survival of bacterial cells within D. 
melanogaster flies after 1, 3 and 5 days of infection. Assays were 
carried out in quadruplicate for two independent experiments and the 
results are expressed as means ± SEM. Statistically significant 
differences at P<0.005 and P<0.05 are identified by ** and *, 
respectively (one-tailed Student's f test). (C) Percentage of flies that 
survived at each time point after being fed with the different P. 
aeruginosa strains. Each point of the survival curve is the average of 
triplicate experiments expressed as means ± SEM. A representative 
experiment out of 3 is shown in the figure. 
doi:1 0.1 371/journal.pgen.1 0041 63.g007 



observed for fldP, whereas transcripts of PA14_22500 showed only 
a small increase (Figure 9B). This result is in line with the 
transcriptional organization of RGP32 described above, and raises 
the question as to how this yZaf-containing operon could be 
regulated. It has been previously observed that genomic islands 
often possess their own regulatory elements. In this sense, gene 
PA14_22550, which is predicted to encode a LysR family 
transcriptional regulator, appears as a promising candidate to 
fulfill this role. We investigated this possibility by using a 
PA14_22550 null mutant strain (ID53714) [28], and compared 
the expression of fldP and its response to H 2 0 2 with those 
previously observed for the isogenic wt strain. Surprisingly, RT- 
PCR analyses showed that inactivation of PA14_22550 produced 
a 3.4-fold increase in the expression of JldP (P= 0.0285), even in 
untreated cells (Figure 6). In fact, treatment with H 2 0 2 only raised 
this difference to 4.1-fold (P— 0.004). Thus, following the increase 
observed due to PA14_22550 inactivation, exposure to H 2 0 2 did 
not produce any further rise of fldP expression (P= 0.4558). Then, 
the collected results indicate that JldP is part of a polycistronic 
operon which is repressed by the LysR-like transcriptional 
regulator present in the same RGP. 

Discussion 

Even under optimal growth conditions a low percentage of the 
electrons involved in cellular redox pathways are diverted to oxygen 
with concomitant ROS generation. This fraction can be dramat- 
ically increased under adverse environmental situations, leading to a 
condition known as oxidative stress. Moreover, cellular auto- 
oxidations are not the only source of ROS and oxidative stress; they 
can also be generated by external redox processes, such as 
phagocytes and other eukaryotic cells, which douse invading 
pathogens with H 2 0 2 as a strategy to prevent infection (reviewed 
in [1,2]). 

A common strategy for coping with ROS damage and oxidative 
stress is the use of electron shuttles to relieve the excess of reducing 
equivalents and redox-active compounds. One of the most 
conspicuous among them is the electron carrier flavoprotein Fid, 
which has been consistently associated with stress protection in a 
number of organisms. The long-chain Fids from Anabaena and E. coli 
(IsiB and FldA, respectively) are encoded by oxidant-inducible 
genes, and confer increased tolerance to environmental and 
nutritional hardships [3,4] . 

Despite the physiological importance of Fid as adaptive resource, 
little had been described about these flavoproteins in the versatile 
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Figure 8. fldP\s a component of the P. aeruginosa variable genome. The fldP gene (in red) is located in a cluster of six genes which is referred 
to as RGP32 [33], here represented in blue in strains PA14, 39016 and HexIT. The genomic context of RGP32 was outlined in gray in the three strains 
and also in the prototypic strain PA01, which does not harbor RGP32. 
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opportunistic human pathogen P. aeruginosa. In this work, we report tire 
presence of a long-chain Fid (FldP) in P. aeruginosa and describe its role 
in the defense of the bacterial cell against oxidative conditions. FldP is 
encoded by the PA14_22540 gene as a 184-amino acid product 
sharing overall sequence similarity and com- 
mon structural signatures with well-characterized flavodoxins such as 
FldA and IsiB (Figure 1). The JldP gene was cloned, expressed in E. coli 
and the resulting product purified to homogeneity, displaying the 
spectral properties and activity of a functional Fid (Figure 2). 

We deemed it necessary to confirm that FldP was a functional 
flavodoxin because short-chain Fids have been already identified 
in both P. aeruginosa and P. put.ida [17,18]. These genes have been 
annotated as mioC by comparison with their ortholog in E. coli. 



Pseudomonas MioC behaves as a functional Fid, mediating FNR- 
catalyzed electron transfer to cytochrome c with a k CM of about 
6 min -1 [18], comparable to that displayed by FldP (22 rnin - , 
Figure 2B), but its physiological role remains yet to be determined. 
Null mutants in mioC did not display growth phenotypes in E. coli 
[37,38] or P. aeruginosa [17]. In the latter organism, however, mioC 
mutations exhibited other pleiotropic effects, including altered 
response to iron stress, increased production of the extracellular 
pigments pyocyanin and pyoverdine, and modified resistance to 
antibiotics [17]. In contrast to the enhanced sensitivity of JldP 
mutants toward H 2 0 2 toxicity (Figure 3A), P. aeruginosa cells 
deficient in mioCwere more tolerant than the wt to oxidative stress 
caused by methyl viologen or H 2 0 2 [17]. 
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Figure 9. Transcriptional organization of RGP32. (A) The fWPgene is co-expressed with genes PA14_22530, PA14_22520 and PA14_22510 as a 
transcriptional unit, thereby constituting an operon (in red). Genes which are not co-transcribed with fldP are shown in blue. The gels below 
represent typical PCR amplifications obtained using cDNA or genomic DNA as templates. (B) Semi-quantitative RT-PCR was performed with specific 
primers for PA14_22530, which shows a H 2 0 2 -induced expression similar to that of fldP. This induction was not observed in PA14_22500, indicating 
that expression of this gene is independent from that of fldP. The rpoD transcripts were used as housekeeping controls. Bars represent relative 
expression levels averaged over two independent experiments (means ± SEM), with expression of the untreated (NT) wt strain normalized to 1. 
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The presence of different non-redundant Fids in a single 
organism is not uncommon. The E. coli genome, for instance, 
contains at least four genes predicted to encode flavodoxins: fldA, 
fldB, mioC, and jqcA [38]. E. coli Fids engage in different cellular 
pathways, and in general they cannot be functionally exchanged 
[38,39]. Moreover, the JldA gene is specifically induced by H 2 0 2 
and redox-cycling oxidants [40,41], whereas mioC is not [37]. 
Likewise, FldP and MioC, while displaying essentially the same 
activity in vitro, appear to play different roles in P. aeruginosa, with 
FldP contributing to the protection against oxidative challenges, 
and MioC in the response to iron stress [17]. 

A search of fldP orthologs in the Pseudomonas genomes 
available in the Pseudomonas Genome Project (http://www. 
pseudomonas.com) [34] and the environmental HexlT strain 
(Feliziani et ai, unpublished) revealed their presence in just six 
out of fourteen strains. The P. aeruginosa genome is made up of 
a conserved core component disrupted by numerous strain- 
specific RGPs, the sum of which conform the accessory and 
variable genome [31,33]. The fldP gene is part of the region of 
genome plasticity RGP32, belonging to the accessory genome 
(Figure 8). Synteny and genome location are highly conserved 
in the PA14, PA39016, NCGM2.S1, B136-33 and HexlT 
strains, with flanking genes being present in PAO 1 , suggesting 
that RGP32 is the result of a common insertion event. In this 
sense, we identified conserved palindromic inverted repeats 
flanking RGP32 in all these genomes, which provide circum- 
stantial evidences on the origin and acquisition of this genome 
block. 



While the conserved repertoire of the core genome codes for 
central functions required for survival and reproduction in any 
habitat [3 1 ,42] , the variable regions may play customized roles in 
adaptation to particular environments [33,43]. This singular 
genome architecture, combining conserved and variable compo- 
nents, bestows P. aeruginosa pangenome upon its capability to 
handle a broad metabolic potential in order to adapt to the widest 
range of environmental niches. In this context, fldP, as part of the 
accessory genome, may provide specialized oxidative stress-related 
functions that benefit survival under stressful conditions, confer- 
ring RGP32 a role as adaptive island. 

Indeed, several lines of evidence indicate that FldP is involved in 
protection against oxidative stress, including (i) the strong 
induction of the fldP gene in response to H 2 0 2 (Figure 6), (ii) the 
enhanced ROS build-up and lower survival of fldP null mutants 
exposed to H 2 0 2 (Figure 3), and (iii) the partial protection 
conferred by FldP overexpression to P. aeruginosa cells deficient in 
mutT against the deleterious effects (Figure 4) and increased 
mutational burden (Figure 5) caused by H 2 0 2 treatment. 

P. aeruginosa houses a multifaceted antioxidant response, 
comprising scavenging enzymes such as catalases and superoxide 
dismutases, thioredoxins and glutaredoxins, as well as small 
antioxidant molecules such as glutathione and melanine [44]. 
OxyR can be considered the master P. aeruginosa oxidative stress 
adaptive response [27]. Furthermore, recent studies have extended 
the implication of OxyR in other P. aeruginosa important responses, 
such as quorum sensing regulation, iron homeostasis and oxidative 
phosphorylation, also identifying a large number of target genes 
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and revealing a far more complex cellular response than previously 
envisaged [45]. It was therefore tempting to speculate that JldP 
could be just another effector gene of this transcriptional regulator. 
Analysis of JldP expression in two oxyR loss-of-function mutants 
ruled out this possibility (Figure 6), indicating that induction oifldP 
in response to oxidative stress proceeds through an OxyR- 
independent pathway. 

Inspection of the gene content and organization of RGP32 
suggested another attractive possibility. The JldP gene is the first in 
a row of five ORFs displaying the same transcriptional orientation 
and extending up to the left end of RGP32 (Figure 8). This 
accessory genomic region contains still another ORF 
(PA14_22550), which is divergently transcribed and encodes a 
putative protein related to the LysR family of transcriptional 
regulators (Figure 8). LysR is the prototype for the most extended 
family of transcription factors in the bacterial world. Although 
originally described as activators of divergently transcribed genes, 
subsequent research placed LysR proteins as global transcriptional 
regulators, acting as either activators or repressors of single or 
operonic genes (reviewed in [46]). We found that four of the five 
divergent ORFs (includingyZrfP) were co-transcribed in response to 
oxidative stress (Figure 9A), therefore constituting an operon. 
Moreover, inactivation of PA14_22550 led to full induction of JldP 
in the absence of oxidants, and to H 2 0 2 insensitivity (Figure 6), 
indicating that the LysR homologue acts as a repressor of fldP and 
presumably the entire operon. The sixth ORF of RGP32, whose 
presumptive product is homologous to protein-disulfide isomerases 
(Table SI), is not part of the operon. 

Stress-associated traits are often encoded by loci adjacent to those of 
other defensive products, especially when they are co-regulated 
[31,42]. In RGP32, the ORF immediately downstream bom JldP 
encodes a putative glutathione S- transferase (Table SI). Involvement 
of this superfamily of conjugative enzymes in the protection against 
oxidative stress has been extensively documented in all types of 
organisms (see, for instance [47]). On the other hand, nucleoside- 
disphosphate-sugar epimerases as that presumably encoded by 
PA14_22520 have been consistently associated with stress responses 
in plants [48], fungi [49] and bacteria [50]. Finally, PA14_22510 
encodes a putative H protein from the glycine cleavage system (Table 
SI). These are lipoate-containing carrier proteins which provide 
reducing power (in the form of thiols) to the multi-enzymatic complex 
involved in glycine decarboxylation. The gene encoding this protein 
has been shown to be strongly up-regulated by H 2 0 2 in streptococci 
[51], and lipoic acid is known to participate in the cellular response 
against oxidative stress in dilferent organisms [52]. Then, all members 
of the operon have the potential to contribute to the defense against 
oxidative stress at different levels. We therefore propose that RGP32 
represents a stress-inducible, self-regulated genetic element which 
confers increased tolerance to oxidative challenges. Under normal 
growth conditions, expression of RGP32 genes should be repressed by 
the LysR-type regulator encoded by PA14_22550. Oxidants 
somehow inactivate this transcription factor allowing induction of 
JldP and other components of the operon. The mechanism by which 
oxidants modulate the activity of the LysR-like protein of RGP32 is at 
present unknown and deserves further investigation. 

The protective effect displayed in vitro by FldP against oxidative 
stress prompted us to evaluate whether this tolerance could be 
advantageous to P. aeruginosa cells exposed to the oxidative assault 
of cells belonging to the mammalian immune system. Indeed, wt 
and complemented P. aeruginosa strains expressing FldP from the 
chromosome or a plasmid exhibited better survival within infected 
macrophages relative to an isogenic mutant lacking this flavopro- 
tein (Figure 7A). Cifani et al. [53] have recently shown that the 
oxidative burst produced by P. aeruginosa-infected macrophages 



plays a key role in the short-term killing of intracellular bacteria 
following invasion, strongly suggesting that the protective effect of 
FldP stems from its antioxidant function as it occurs in vitro. 

We also used the insect model system of D. melanogaster (in which 
the PAH strains has been shown to be particularly aggressive 
[29]), to further investigate the importance of this protective effect 
on a real infection process. In good agreement with the differential 
sensitivity observed in macrophages, JldP mutant bacteria accu- 
mulated to lower levels in D. melanogaster, as compared to the wt or 
the complemented strains (Figure 7B). Noteworthy, the mortality 
kinetics was delayed for ~24 h in the mutant-infected flies, 
whereas the death rate of the flies infected with the yWP-deficient 
bacteria complemented with p2-JldP was equivalent to that of the 
wt strain (Figure 7C). Thus, the oxidative killing of P. aeruginosa 
within Drosophila hemolymph may involve mechanisms similar to 
those utilized by mammalian hosts. While the FldP effect suggests 
that the electron shuttle could be involved in P. aeruginosa virulence, 
it is more likely that the different death rates observed in Figure 7C 
are a consequence of longer persistence of FldP-containing 
bacteria in the fly due to the adaptive advantages conferred by 
the flavoprotein. In line with this proposal, mutation of the P. 
aeruginosa oxyR gene had similar effects on bacterial survival and 
Drosophila killing as those reported here [54] . 

Then, our data identify an oxidant-responsive long-chain 
flavodoxin in P. aeruginosa, which participates in the defense against 
ROS and contributes to the bacterial tolerance to the oxidative 
onslaught elicited by the host immune system. We anticipate that 
studies on this direction will lead to a more comprehensive 
panorama of the mechanisms allowing this opportunistic pathogen 
to adapt and persist in stressful and dynamic environments. 

Materials and Methods 

Bacterial strains, plasmids and growth conditions 

P. aeruginosa PA 14 and its isogenic strains ID38939 
(PA14_22540 mutant, JldP), ID53714 (PA14_22550 mutant) and 
ID54029 (PA14_70560 mutant, oxyR) were kindly provided by Dr 
Eliana Drenkard and Dr Jonathan Urbach from the Massachusetts 
General Hospital, Boston, USA [28], whereas P. aeruginosa 
MPAOl and its isogenic mutT strain were provided by Dr 
Michael Jacobs from the University of Washington Genome 
Center, USA [55]. Insertion of the MAR2xT7 mini-transposon in 
mutant ID38939 is unlikely to have polar effects on the expression 
of downstream genes since the consecutive aacCl promoter is 
oriented in the same direction. To prepare inocula, bacteria were 
routinely cultured on Luria-Bertani (LB) agar plates from frozen 
stocks and subcultured overnight in LB liquid medium at 37°C 
with shaking at 220 r.p.m. Antibiotics were used at the following 
concentrations: 30 u,g ml -1 gentamicin (Gm); 250 |J.g ml -1 kana- 
mycin (Km). 

Bioinformatics analysis 

To find Fid homologs in P. aeruginosa PA 14, the Anabaena IsiB 
and E. coli FldA sequences were compared against the entire PA 14 
genome using the DELTA-BLAST Search Tool [19]. Multiple 
DNA sequence alignments were performed by using ClustalW 
(http://www.clustal.org). Secondary structures were predicted 
using the Jpred3 software provided by the Dundee Scotland 
University (http://www.compbio.dundee.ac.uk/ www-jpred). 

Cloning of the fldP gene from P. aeruginosa PA14 and the 
isiB gene from Anabaena PCC7119 

A DNA fragment containing the entire coding region of the JldP 
gene (PA14_22540) was amplified by PGR from PA14 genomic 
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DNA using oligonucleotides FldP-F and FldP-R (Table S2), 
containing BamHI and Hindlll restriction sites, respectively. The 
PGR product was ligated to the pGem-T Easy vector (Promega) 
and subsequently cloned into the broad-host-range plasmid 
pBBRlMCS2 (p2), which harbors a Km resistance marker [56], 
to generate p2-fldP. A similar strategy was employed to prepare 
p2-isiB containing the Fid-encoding gene from Anabaena PCC7 1 19, 
originally cloned in pEMBL8-wtB [20], using oligonucleotides 
IsiB-F and IsiB-R as forward and reverse primers, respectively 
(Table S2). The resulting plasmids (p2-jldP and p2-isiB) and the 
empty p2 vector (as control), were introduced in the different P. 
aeruginosa strains via electroporation [57]. 

Construction of an oxyR deletion-mutant strain 

Complete deletion of the oxyR gene was carried out as previously 
described [57]. All primer sequences are described in Table S2. 
Briefly, a first round of three PCR reactions was performed in 
which the 5' and ?>' flanking regions of oxyR, as well as a Gm 
resistance cassette were amplified from plasmid pPS856 [58] using 
four gene-specific primers (Oxy-UpF-GWL, Oxy-UpR-Gm, Oxy- 
DnF-Gm and Oxy-DnR-GWR) and the common Gm-specific 
primers (Gm-F and Gm-R). This generated three fragments with 
partial overlaps either to each other or the attBl and attB2 
recombination sites. The purified fragments were then assembled 
in vitro by overlap extension during the second round PCR using 
the common primers GW-attBl and GW-attB2. This resulted in 
an ory-ff-deletion-mutant PCR fragment which was subsequendy 
cloned into pDONR221 (Invitrogen) via the BP clonase reaction 
to create pDONR221-oxiJ-fl::Gm. This construct served as the 
substrate for LR clonase-mediated recombination into the 
destination vector pEX18ApGW. The resulting suicide vector 
pEXISApGW-oxy-fciGm was then transferred to P. aeruginosa and 
the plasmid-borne oxjiff-deletion mutation was exchanged with the 
chromosome via homologous recombination to generate the 
chromosomal deletion mutant. 

Expression and purification of recombinant FldP and IsiB 

For expression in E. coli, a 568-bp fragment encoding the 
complete sequence of the JldP gene was obtained by PCR 
amplification, using p2-JldP as template, and primers Rec-FldP-F 
and Rec-FldP-R, which contain restriction sites for Mdel and 
Hindlll, respectively (Table S2). The amplified fragment was 
digested with the corresponding enzymes, cloned into compatible 
sites of pET-TEV (Novagen) under the control of the T7 
promoter, and fused in-frame to an N-terminal His-tag. To 
improve solubility, BL2 1 E. coli cells were co-transformed with this 
vector and plasmid pG-Tf2 (Takara Bio Inc) expressing E. coli 
molecular chaperones (GroEL, GroES and Trigger Factor). After 
induction with 0.2 mM isopropyl-fi-D-thiogalactoside (IPTG), the 
soluble flavoprotein was purified from cleared lysates in a Ni-NTA 
column by elution with 500 mM imidazole. Expression and 
purification of IsiB were carried out according to Fillat et al. [20]. 
UV-visible spectra of the purified recombinant proteins and FMN 
were recorded in 50 mM Tris-HCl pH 8.5. 

Functional characterization of FldP 

The ability of both electron shuttles, FldP and IsiB, to mediate 
the cytochrome c reductase activity of Anabaena FNR was assayed 
according to Shin [59]. The reaction mixture contained 3 mM 
glucose 6-phosphate, 0.3 mM NADP + and 1 unit ml 1 glucose 6- 
phosphate dehydrogenase (G6PDH, to generate NADPH), 
0.5 |J.M FNR, 50 |jM equine heart cytochrome c and various 
amounts of Fid in 50 mM Tris-HCl pH 8.5. Cytochrome c 



reduction was followed at 30°C by the increase in absorbance at 
550 nm (e 55 o= 19 mM -1 cm -1 ). 

ROS detection 

A modified version of the FOX II assay [60] was used to 
quantify the presence of peroxides in bacterial extracts. Cultures of 
the parental PA 14 and the JldP mutant strains transformed with 
either p2 or p2-fldP were grown aerobically in LB broth at 37°C 
for 5 h with the appropriate antibiotics. Then, H 2 0 2 was added to 
final concentrations of 0, 25 and 50 mM, and bacterial 
suspensions were incubated for 30 min with vigorous shaking. 
Cultures were split into two equal portions; one of them was used 
to measure protein concentration, while the other was centrifuged, 
washed with 0.9% (w/v) NaCl and finally resuspended in 1 ml of 
an 80:20 ethanol/water solution containing 0.01 % (w/v) butylated 
hydroxytoluene (BHT). Samples were disrupted by sonic oscilla- 
tion (10 times for 10 sec each, 30% amplitude), centrifuged at 
10,000 g for 10 min, and 250 |J,1 of the supernatants were 
combined with 250 ju.1 of 10 mM Tris-phenyl phosphine in 
methanol (TPP, a -OOH reducing agent), or with 250 Jul of 
methanol, to measure total oxidants. Mixtures were incubated for 
30 min to allow complete -OOH reduction by TPP. Five hundred 
ul of FOX reagent (100 uM xylenol orange, 4 mM BHT, 250 mM 
ferrous ammonium sulphate and 25 mM H 2 S0 4 in 90% (v/v) 
methanol) were then added to each sample, and the absorbance at 
560 nm was recorded 10 min after reagent addition. The 
absorbance differences between equivalent samples with and 
without TPP indicate the amounts of -OOH, which were 
calculated using a 0-20 uM H 2 0 2 standard curve. Protein 
concentrations were estimated in cleared lysates in 50 mM Tris- 
HCl pH 8.0 by a dye binding method [61], using bovine serum 
albumin as standard. 

For observation at the confocal microscope, DCFDA was 
introduced into P. aeruginosa cells by electroporation. Briefly, 10-ml 
overnight cultures of the various strains were collected by 
centrifugation at 12,000 g for 2 min, washed three times with 
1 ml of 0.3 M sucrose and finally resuspended in 100 ul of the 
same solution. DCFDA (in dimethyl sulfoxide) was added to a final 
concentration of 500 uM and the suspension transferred to 
electroporation cuvettes of 1-mm width. Cells were electroporated 
in an Electro Cell Manipulator 600, BTX electroporation System 
at 25 uF, 2.5 kV and 200 Q. for 5 ms, diluted with 300 ul of 0.3 M 
sucrose and divided into two equal fractions. One of them was 
incubated with 25 mM H 2 0 2 at 37°C for 15 min and the other 
kept under the same conditions without the oxidant. Seven |il of 
the suspensions were mixed with 3 |J.l of the membrane-staining 
dye FM 4-64 (0.1 mg ml ) on the surface of a glass plate and the 
lid was sealed with colorless nail paint. Cells were observed at 
excitation/emission maxima of ~5 15/640 nm in a Nikon TE- 
2000-E2 confocal microscope. The number of fluorescent spots 
was determined for each dye by using the Image J 1.46 software. 

H 2 0 2 susceptibility measurements 

For H 2 0 2 susceptibility assays, cells were grown to mid- 
exponential phase (OD 600 ~0.8), collected by centrifugation and 
washed with 0.9% (w/v) NaCl solution. Then, bacteria were 
normalized to ~10 cells and treated with 50 mM H 2 0 2 for 
30 min at 37°C. Cells were centrifuged again, washed twice with 
0.9% (w/v) NaCl and finally resuspended in fresh LB medium. 
Serial dilutions were spotted on LB agar plates and incubated 
overnight at 37°C to determine viability. Non-treated controls 
were included. Bacterial susceptibility was expressed as the ratio 
between the survivals of treated to untreated cells. 
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Estimation of spontaneous and H 2 0 2 -inclucecl mutant 
frequencies 

For estimation of spontaneous-mutation frequencies, five bacterial 
colonies of the different P. aeruginosa strains were cultured overnight in 
10 ml LB medium at 37°C with shaking at 220 r.p.m. Appropriate 
dilutions of the cultures were plated on LB agar to determine the total 
number of viable cells, or on LB agar supplemented with 
500 ug ml 1 streptomycin to count the number of streptomycin- 
resistant cells, after overnight incubation at 37°C. Spontaneous- 
mutation frequency was determined as the ratio between the number 
of streptomycin-resistant cells and the number of viable cells. 

Determination of H 2 O z -induced mutant frequency was carried 
out according to previously described protocols [62], with some 
modifications. Briefly, five independent bacterial cultures for each 
strain were grown in LB to mid-exponential phase, collected by 
centrifugation and washed with 0.9% (w/v) NaCl. Cells were 
normalized to ~10 8 and subsequently treated with 50 mM H 2 O z 
for 30 min at 37°C, centrifuged again and washed twice with 1 ml 
of 0.9% (w/v) NaCl. An aliquot of treated cells (0.5 ml) was diluted 
10-fold in fresh LB broth and cultured overnight at 37°C with 
shaking at 220 r.p.m. Then, 100 ju.1 of each overnight culture were 
plated onto LB agar supplemented with 500 u,g ml 1 streptomy- 
cin, whereas aliquots from appropriate dilutions were plated on LB 
agar without antibiotic to measure cell viability. The H 2 0 2 - 
induced mutation frequency was calculated as above. 

Semi-quantitative RT-PCR 

P. aeruginosa PA14 mid-exponential phase cultures (OD 600 ~0.8) 
were treated with 50 mM H 2 0 2 for 15 min, centrifuged and 
subsequently washed twice with 0.9% (w/v) NaCl. Treated and 
untreated cultures were used to extract total RNA using the RNA 
Purification Kit (Fermentas). RNA was quantified by UV spectro- 
photometry, and its integrity was checked by electrophoresis in 
1.5% (w/v) agarose gels. Then, 1 ug of total RNA was reverse- 
transcribed using the QuantiTect Reverse Transcription Kit 
(QIAGEN). PCR primers were manually designed with the 
assistance of the Netprimer software (PREMIER Biosoft Interna- 
tional, Palo Alto, CA) and evaluated for their specificity with the 
BLAST program at the NCBI Web site. Specific transcripts were 
semi-quantitatively measured by RT-PCR using primers FldP- 
RT-F and FldP-RT-R (for JldP), 22530-RT-F and 22530-RT-R 
(for PA14_22530), and 22500-RT-F and 22500-RT-R (for PA14_ 
22500). Transcripts of the rpoD gene were amplified with primers 
RpoD-RT-F and RpoD-RT-R and served as housekeeping con- 
trols. All primer sequences are described in Table S2. The optimal 
number of cycles was determined in advance to evaluate expres- 
sion in the exponential phase of amplification. Final cycling 
conditions included a hot start at 95°C for 10 min, followed by 32 
cycles at 94°C for 30 sec, 60°C for 30 sec and 72°C for 30 sec, 
and a final extension at 72°C for 5 min. Specificity was verified by 
agarose gel electrophoresis. Fold change in gene expression was 
calculated by measuring band intensities with the Gel-Pro An- 
alyzer Software. No amplification was observed in PCR reactions 
containing water or non-reverse transcribed RNA as template. 

Characterization of RGP32 transcriptional organization 

Primers were designed in order to amplify PCR products 
containing regions of two neighboring genes for those genes of 
RGP32 which share the same orientation and are located 
downstream from JldP (Figure 9A). The cDNA to be used as a 
template for PCR was obtained by reverse transcription of purified 
total RNA as described above. Thus, primers FldP-F and FldP-30- 
R were employed to determine the co-expression of genes JldP and 



PA14_22530 (Table S2). In the same way, co-expression of 
downstream contiguous genes was analyzed with the following 
primers: a) 30-20-F and 30-20-R to analyze PA14_22530 and 
PA14_22520 co-expression; b) 20-10-F and 20-10-R for 
PA14_22520 and PA14_22510; and c) 10-00-F and 10-00-R for 
PA14_22510 and PA14_22500 (Table S2). Genomic DNA was 
used as template for positive controls of the PCR reactions. No 
amplification was observed in PCR reactions containing water or 
non-reverse transcribed RNA as template. 

Intracelullar survival in macrophagic cells 

Cells of the macrophagic-derived line RAW 264.7 were grown on 
1 2-mm-diameter wells until they reached 95-100% confluence. To 
prepare inocula, overnight cultures of P. aeruginosa were washed with 
phosphate buffered saline (PBS) and suspended in Dulbecco's 
Modified Eagle's Medium (DMEM) at a final concentration of 10 7 
cells ml 1 . Bacterial antibiotic protection assays were conducted as 
previously described with few modifications [63]. Briefly, macro- 
phagic cells were inoculated with the various P. aeruginosa strains by 
the addition of 500 ui of the bacterial suspension, corresponding to 
a multiplicity of infection (MOI) of about 20, followed by 
centrifugation for 10 min at 1000 r.p.m. to facilitate cell contact. 
After a 2-h incubation period at 37°C, the supernatant was removed 
and the cells were washed three times with PBS to remove non- 
associated bacteria. To count intracellular P. aeruginosa, fresh 
DMEM with streptomycin (1 mg ml -1 ) and carbenicillin 
(400 |J.g ml" 1 ) was added to the cell monolayer and incubated for 
2 h to kill extracellular bacteria. Antibiotic bactericidal activity was 
confirmed by plating 1 00-ixl aliquots from the wells directly on LB- 
agar. Cells were again washed to remove antibiotics, then lysed by 
incubation with 0. 1 % (v/v) Triton X- 1 00 in PBS. At this initial stage 
(Ti) intracellular bacteria were enumerated by plating serial 
dilutions of cell lysates on LB agar and counting CFU. To estimate 
P. aeruginosa intracellular survival, an equivalent set of wells were 
incubated with antibiotics for an additional 3-h period (Tfj, cells 
were lysed and surviving bacteria enumerated as described above. 
Intracellular survival was calculated as the ratio between the final 
and the initial (Tf/Ti) CFU counts. Three wells of cells were used 
for each strain in each experiment, and all experiments were 
repeated three times. To rule out the pre-existence of resistant 
bacteria in the overnight cultures, 1 0 x inocula (5 x 1 0 7 cells) were 
plated on LB agar containing streptomycin (1 mg ml ) and 
carbenicillin (400 |J,g ml" ), with not a single CFU observed even 
after 72 h of incubation. 

Phagocytes viability was monitored at every time point of the 
experiment by measuring release of LDH (CellTiter 96 Aqueous 
Non-Radioactive Cell Proliferation Assay, Promega), which 
showed an initial decrease of ~50% in the first 2 h, but remained 
unaltered after the antibiotic treatment. 

Bacterial survival in Drosophila melanogaster 

P. aeruginosa cells of the various strains were grown overnight, 
washed and diluted in PBS, 10% (w/v) sucrose, 250 ug ml" 1 Km 
to OD 6( ,o = 0.025. D. melanogaster strain Wl 1 18 flies (4-5 days old) 
were starved for 3 h and then fed continuously at 25°C on cottons 
plugs, which had been previously embedded in the bacterial 
solution. For each time point, 3 groups of 5 flies were pestle 
homogenized in 200 ju.1 of PBS and the homogenate was serially 
diluted in PBS and plated on LB agar containing 
250 |Xg ml" 1 Km, to quantify the number of bacterial colonies. 

Fly survival experiments 

Overnight cultures of each P. aeruginosa strain were diluted in 
fresh LB to an OD 60 o = 0.25. This solution was then diluted 10- 
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fold in PBS, 10% (w/v) sucrose, 250 u.g ml 1 Km. Groups of 15 
D. melanogaster Wl 1 18 male flies were starved during 3 h and then 
placed in vials with sterile cotton plugs, which had been previously 
embedded in 3 ml of the bacterial suspension. Flies were kept at 
25°C and survival was monitored daily. Three groups of 15 flies 
were used for each condition in three independent experiments. 

Statistical analysis 

Statistical analyses were performed using one-tailed Mann- 
Whitney test appropriate for nonparametric adjustment. When 
appropriate, one- or two-tailed Student's t test was applied. In all 
cases, P values less than or equal to 0.05 were considered 
statistically significant. 

Ethics statements 

The collection of P. aeruginosa clinical isolates used in this work 
has already been published by Feliziani et al. [36]. In this previous 
study, the informed consent as well as the approval from an 
Institutional Research Committee were appropriately evaluated 
and fulfilled the PLOS ethical standards. 

Supporting Information 

Figure SI Expression and purification of FldP from E. coli 
extracts. (A) FldP was recovered as insoluble inclusion bodies 
even after expression at low temperature and IPTG concentra- 
tion. BL21 E. coli cells transformed with plasmid PET-TEV 
carrying the jldP gene were grown in LB medium at 20 or 37°C 
for 3 hours with 10 or 50 |J,M IPTG. Bacteria were incubated at 
the indicated temperature for 2 hours and cells were harvested, 
lysed and centrifuged. Cleared extracts were resolved by SDS- 
PAGE and stained with Coomassie Brilliant Blue. Lanes 1—5, 
supernatants corresponding to 20 p.g of soluble protein; lanes 6- 
10, pellets representing an equivalent amount of cells. Lanes 1 
and 6, without IPTG; lanes 2-3 and 7-8 with 10 uM IPTG at 20 
and 37°C, respectively; lanes 4-5 and 9-10 with 50 uM IPTG at 
20 and 37°C, respectively. (B) Purification of FldP. After 
induction of expression with 0.2 mM IPTG in E. coli BL21 
expressing the chaperones GroEL, GroES and Trigger Factor, a 
minor soluble fraction of the flavoprotein was isolated by passage 
through a Ni-NTA agarose column (QIAGEN). Purified FldP 
migrated as a polypeptide of ~24 kDa. Lanes 1-2, pellet and 
supernatant, respectively corresponding to 20 ug of soluble 
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(TIF) 
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